Extraction of domain-specific magnetization reversal for nanofabricated periodic 
arrays using soft x-ray resonant magnetic scattering 
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A simple scheme to extract the magnetization reversals of characteristic magnetic domains on 
nanofabricated periodic arrays from soft x-ray resonant magnetic scattering (SXRMS) data is pre- 
sented. The SXRMS peak intensities from a permalloy square ring array were measured with field 
cycling using circularly polarized soft x-rays at the Ni L3 absorption edge. Various SXRMS hys- 
teresis loops observed at different diffraction orders enabled the determination of the magnetization 
reversal of each magnetic domain using a simple linear algebra. The extracted domain-specific hys- 
teresis loops reveal that the magnetization of the domain parallel to the field is strongly pinned, 
while that of the perpendicular domain rotates continuously. 

PACS numbers: 75.25.+Z, 75.75.+a, 75.60.-d 
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Understanding the reversal mechanism of the magne- 
tization in periodic arrays of submicron and nanoscale 
magnets is of both scientific and technological interest. 
Fundamental changes in the statics and the dynamics 
of magnetization reversal imposed by nanostructures en- 
rich the physics of nanomagnetism. A precise control 
of magnetization reversal involving well-defined and re- 
producible magnetic domain states in nanomagnet ar- 
rays is key to future applications, such as high density 
magnetic recording* or magnetoelectronio* devices. To 
achieve this, topologically various nanomagnets, ranging 
from simple disks* to more complicated rings* or negative 
dots (holes)*, have been investigated. However, as rather 
well-defined but non-single magnetic domains (or domain 
states) form in complicated geometries, it becomes diffi- 
cult to characterize precisely magnetization reversal in- 
volving each domain at small-length scales with cither 
conventional magnetization loop measurements, such as 
magneto-optical Kerr effect (MOKE) magnetometry, or 
magnetic microscopy, such as magnetic force microscopy 
(MFM). Moreover, in large-area arrays typically cover- 
ing areas of a few square millimeters, extracting over- 
all domain structures during reversal from microscopic 
images is clearly unreliable. Though diffracted MOKE 
measurement has been proposed recently to deal with 
this problem, it is found to provide little quantitative 
information on magnetization reversal involving domain 
formation and is limited to micrometer-length scales^ 

Such quantitative information is available using the 
technique of soft x-ray resonant magnetic scattering 
(SXRMS). 7 This technique exploits strong enhancement 
of the magnetic sensitivity of scattering intensities when 
incident circularly polarized soft x-rays are tuned to an 
absorption edge of constituent magnetic atoms. SXRMS 
has been used to study the magnetic structure in mag- 
netic thin films* or periodic arrays of stripe domains or 



nanolincs 9 . In this paper we present a simple scheme 
to extract quantitatively domain-specific magnetization 
reversal for nanomagnet arrays from SXRMS measure- 
ments. For this purpose, sample rocking curves, yielding 
in-plane diffraction scans, have been measured and an- 
alyzed on the basis of our previous work performed in 
the hard x-rays^* In order to obtain magnetic informa- 
tion, SXRMS peak intensities were measured by vary- 
ing the applied field at different diffraction orders, whose 
scattering structure factors are different. This allows us 
to determine directly the magnetization reversal of each 
magnetic domain using a simple linear algebra. The basic 
idea of incorporating such nonuniform magnetic domains 
in scattering theory has been explored in our previous 
work on polarized neutron scattering^* 

For this study an array of permalloy (Ni§oFe2o) square 
rings was fabricated by a combination of e-beam lithogra- 
phy and lift-off techniques. A standard silicon wafer was 
spin-coated with a double-layer positive-type e-beam re- 
sist, and the resist layer was then patterned by e-beam 
lithography. A 20-nm-thick permalloy film was deposited 
onto it using an electron-beam evaporator in a vacuum of 
about 10 -8 Torr. The as-deposited unpatterned film was 
magnetically soft with coercive and uniaxial anisotropy 
fields of a few Oersteds. Finally, after ultrasonic-assisted 
lift-off, the square rings were arranged in an array of 2 x 
2 mm 2 . 

Experiments were performed at sector 4 of the Ad- 
vanced Photon Sourceii* Polarized soft x-rays at the 
beamline 4-ID-C were generated by a novel circularly po- 
larized undulator that provided left- and right-circular 
polarization switchable on demand at a polarization > 
96%. The photon energy was tuned to the Ni L3 absorp- 
tion edge (853 eV) to enhance the magnetic sensitivity. 
While a vacuum compatible sample stage was rotated, 
the diffracted soft x-ray intensities were collected by a 
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FIG. 1: Diffraction intensities of the rocking scan along the 
q x direction at q z = 0.955 nm -1 from a square ring array. 
Circularly polarized soft x-rays were used and the photon en- 
ergy was tuned to the Ni L3 absorption edge (853 eV). Circles 
represent measurements, and lines represent the calculations. 
Insets: (a) scanning electron micrograph and (b) schematic 
of the square ring array studied. The scale bar is 1 /im. (c) 
and (d) show schematics of the scattering geometry, where 6i 
and Of are the incidence and exit angles, and H and M are 
the applied field and magnetization, respectively. 



Si photodiode detector with a fixed angle of 12.8°. This 
sample rocking scan, where the incident and exit angles 9i 
and 8f were varied with the total scattering angle (0i+0f) 
fixed, yielded a transverse q x scan at a fixed q z value (see 
Fig. 1). The angular resolution was defined by a pinhole 
between the sample and detector to be about 0.03°. The 
sample was mounted in the gap of an eletromagnet that 
provides fields in the scattering plane of up to ±800 Oe. 

Figure 1 shows the diffraction intensities of the sample 
rocking scan measured as a function of q x at q z — 0.955 
nnf 1 from the square ring array with the saturation 
field. Diffracted intensities show peaks corresponding to 
an array period of 1.151 /im. Following Refs. [TolfTlL the 
diffracted intensity I can be expressed in the kinematical 
approximation as 

I(q x ;H) = \F(q z )\ 2 £ \ Pc F c (n x ,\n y \) 

Tlx 1 1 fly I 

2 

+ pM F M(n x ,\n v \;H) K(q x ;n x ,\n y \) , (1) 

where H is the applied field, F(q z ) is the ID form factor- 
along q z direction and consequently a constant value for 
a fixed q z , and i*c7(M) an d Pc(M) are the charge (mag- 
netic) form factors on the q x — q y plane and the charge 



(magnetic) contributions to the total atomic scattering 
amplitude, respectively. Near resonance energies pu is 
proportional to the vector product (e^ x e^) -M^ where §j 
and &f are unit photon polarization vectors for incident 
and scattered waves, respectively, and M is the magneti- 
zation vector. For circularly polarized beams, this vector 
product reduces approximately to the component M x in 
the inset (c) of Fig. 1. Therefore, the magnetization 
referred to hereafter represents strictly the parallel com- 
ponent to the x-axis or the field direction (H \\ x in this 
study) of the magnetization vector, indices for 

Bragg points in the reciprocal space with the relationship 
of q XiV = (2Tr/a)n Xt y, where a is the period of the array. 
Since the resolution function 1Z is a long thin ellipse ori- 
ented in the q y direction, nonzero n y values should also be 
taken into account for a g^-scan performed at q y = OiSiiL 
We also note that the evolution of the peak widths of the 
different diffracted orders as a function of q x , as shown 
in Fig. 1, was calculated using the model proposed by 
Gibaud et alM 

For a square ring, the charge form factor Fc(n x , \n y \) 
can be expressed by 



Fc{n x , \n y \) = C 



sine (nxfri + 72 )) sine (jn,, I (71 +72)) 
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, (2) 



where C is 1 for \n y \ — and 2 for \n y \ ^ 0, and 
sinc(a;) = sin(x)/x. 71,2 = 2^1,2/0, where di and di 
are the width of the ring and half of the inner square 
size, respectively, as shown in inset (b) of Fig. 1. Since 
Fq and Fm are the same functional forms for a saturated 
uniform magnetization, the diffracted intensities can be 
calculated using Eqs. (1) and (2) and are shown as the 
solid line in Fig. 1. From the best fit, d\ and di were 
estimated to be 162 ±4 and 377 ±4 nm, respectively, and 
subsequently the gap between rings was 73 nm. 

Figure 2 shows SXMRS peak intensities measured at 
various diffraction orders while field cycling. All field 
dependencies show magnetic hysteresis loops but with 
different features. This is due to different magnetic 
form factors for different diffraction orders, which re- 
flect nonuniform domain formation during magnetization 
reversal, as pointed out in diffracted MOKE studies. 6 
Therefore, the magnetic form factor Fm in Eq. (1) should 
be expressed by the sum of the contributions of all pos- 
sible magnetic domains, i.e., 



F M (n x , \n y \;H) = ^mi{H)F M > (n Xl \n v \), 



(3) 



where mi(H) represents the field-dependent magnetiza- 
tion of the l-th domain, which is the quantity of interest. 
It is noticeable that the magnetic factor in Eq. (3) can be 
factorized into the field- and structure-dependent factors. 
However, the domain-specific magnetization mi(H) can- 
not be directly extracted from measured SXRMS hystere- 
sis loops because, as described in Eq. (1), the diffracted 
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FIG. 2: SXRMS magnetic hysteresis loops (circles) measured 
at several diffraction peaks, indicated by numbers in Fig. 1. 
The solid lines represent the calculated hysteresis loops. 
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FIG. 3: (a) and (d) SXRMS hysteresis loops with positive- 
to-negative intensities (I P n) nipped with respect to H = for 
zeroth- and third-order peaks, respectively. Difference [(b) 
and (e)] and sum [(c) and (f)] intensities between I np (H) 
(negative-to-positive) and nipped I pn (—H) intensities as a 
function of the applied field. The solid lines in (c) and (f) 
represent the calculations. 



intensities are the absolute square of the sum of the struc- 
tural and magnetic contributions. 

To tackle this problem, we considered the difference 
between the field-dependent intensities of I np (H) and 
Ipn{—H), where I np , I pn represent the intensities mea- 
sured while the field is swept along the negative-to- 
positive and positive-to-negative directions, respectively, 
and I pn (—H) represents the intensities flipped from 



Ipn(H) with respect to H = 0, as shown in Fig. 3. 
Assuming that the magnetization reversal has inversion 
symmetry about the origin, only the difference between 
Inp(H) and I pn (—H) is the opposite sign of mi(H) in 
Eq. (3). However, this does not mean that they are 
symmetrical to a certain horizontal line in Fig. 3 (a) 
or (d) because the intensities contain quadratic terms to 
the magnetization mi(H), as described in Eq. (1), con- 
sequently leading to a nonsymmetry to the origin (or the 
center of mass) of SXRMS hysteresis loops in Fig. 2. 
This effect is clearly seen in the sum intensities of Figs. 
3 (c) and (f) and will be an important characteristic of 
scattering-based hysteresis loops. However, in turn, these 
quadratic terms can be ruled out by taking the difference 
intensities, which result, linearly proportional to 

the magnetization. 

These difference intensities at the n x -order peak can 
be expressed from Eqs. (1) and (3) as 



A„ x (if) = I np (n x ;H) - I pn (n x ; 

= me[p c p* M }n n jF(q z 



-H) 



(4) 



Tlx-; fay |) ft* M {P'&i 

I |n„|=0 

where 72.„ x and lZ n represent resolution functions evolv- 
ing n x and n y indices, respectively, into which lZ(n x , \n y \) 
in Eq. (1) can be factorized. If we further normalize A„ x 
by its maximum intensity with a saturation magnetiza- 
tion m s , we can obtain a set of linear equations as 



where 



Br, 
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mi{H) 



ly\)K n 



(5) 



(6) 



Here we used the relationship of Ya* 1 ( n »> \ n v\) = 
Fc\Ji>xi \ n y\) ■ Applying linear algebra, the normalized 
magnetizations mi—x t ...^{H)/m s of N domains can be 
finally obtained directly from the normalized difference 
intensities A„ x (7?)/|A™ ax | measured at N different n x 
orders by taking the inverse of N x N matrix B nic i. 

In principle, the field-dependent intensities measured 
at semi-infinite numbers of orders can thus be used to de- 
termine the magnetization reversal of each infinitesimal 
cell in a unit nanomagnet. However, this is practically 
restricted due to finite measurable diffraction peaks and 
a high symmetry of experimental geometry. The latter 

gives rise to a strong dependence between {n x , \n y \) 
with different domain I or diffraction order n x , and, as 
a result, makes the matrix B n ^i singular and noninvert- 
ible. To lower the geometrical symmetry of experimental 
configuration, we can use a vector magnetometry setup 
by rotating the sample/electromagnet assembly with re- 
spect to the incident photon direction^ Since SXRMS 
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FIG. 4: Extracted magnetization reversals along the field di- 
rection of four types of characteristic domains, which are de- 
picted as gray-filled regions in the insets, from SXRMS hys- 
teresis loops at zeroth-, first-, third-, and fifth-order peaks. 

intensities are proportional to the component of the mag- 
netization vector along the projected incident photon di- 
rection onto the sample surface, this approach can also 
provide vectorial information about domain-specific mag- 
netization. 

In our setup, where both incident beam and field direc- 
tions are parallel to one of the sides of the square rings, 
there may be four characteristic domains, as shown in 
the insets of Fig. 4. Each domain consists of two or 
four subdomains, whose structure-dependent form fac- 
tors in Eq. (3) are identical, and, therefore, its mag- 
netization represents an average value over subdomains. 
The explicit expressions of the structure-dependent mag- 
netic form factors Fj$ of these four domains for (n x , \n y \) 
diffraction order in Eq. (3) will be presented elsewhere^ 
We note that these domains have been chosen to mini- 
mize sigularity of the matrix B rix i and may not be thus 
energetically viable. Nevertheless, this scheme can unam- 
biguously provide information about domain formation 
by considering snapshots of the resultant domain-specific 
magnetizations at each field. 

To construct a 4 x 4 matrix B rix ; for four magnetic 
domains, at least four different peak intensities are re- 



quired, and fifteen combinations can be allowed for six 
measured peaks, as shown in Fig. 2. However, we ex- 
cluded the second peak due to its bad statistics and also 
some other combinations due to relatively small values 
of the determinants of their matrices -B^;, leading to a 
singularity of the matrix. For the optimum combination 
for non-singularity, zeroth-, first-, third-, and fifth-order 
peaks were then chosen. Figure 4 shows the finally ex- 
tracted magnetization reversals for each domain using 
Eqs. (5) and (6). All magnetizations in Fig. 4 are nor- 
malized by the saturation and represent the components 
projected along the field direction, as discussed above, of 
the magnetization vectors. To confirm these results, the 
sum intensities in Figs. 3(c) and 3(f) and the SXMRS 
hysteresis loops for all observed peaks in Fig. 2 were 
also generated by substituting the results in Fig. 4 into 
Eq. (1). These calculations (solid lines) show a good 
agreement with the measurements. 

A remarkable feature in the extracted magnetization 
reversals is that while the domain (I) perpendicular to 
the field rotates coherently, the parallel domain (II) is 
strongly pinned. Interestingly, this is similar to the do- 
main behaviors in the antidot arrays^i^ whose geometry 
resembles the square ring array except for narrow gaps 
between rings. On the other hand, the domain (II) clearly 
shows plateaus, which have been observed generally in 
circular or octagonal ring magnets and are attributed to 
the vortex stated A detailed discussion is beyond the 
scope of this paper and will be presented elsewhere^ 

In summary, we successfully demonstrated that 
domain-specific magnetization reversals can be extracted 
directly from SXRMS hysteresis loops measured at var- 
ious diffraction orders. Extracted domain-specific mag- 
netization reversals are expected to provide a new in- 
sight into magnetic switching mechanism on nanofab- 
ricated arrays. Future studies, exploiting the element- 
selectivity and vector magnetometry setup, will provide 
further three-dimensional information in nanofabricated 
multilayers such as giant magnetoresistance and pseu- 
dospin valve structures. 
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